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Abstract 

In this paper, we develop approximative two-flavor neutrino oscillation formulas including sub- 
leading nonstandard interaction effects. Especially, the limit when the small mass-squared differ- 
ence approaches zero is investigated. The approximate formulas are also tested against numerical 
simulations in order to determine their accuracy and they will probably be most useful in the GeV 
energy region, which is the energy region where most upcoming neutrino oscillation experiments 
will be operating. Naturally, it is important to have analytical formulas in order to interpret the 
physics behind the degeneracies between standard and nonstandard parameters. 
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I. INTRODUCTION 



Neutrino oscillation physics has had a remarkable progress as the leading description for 
neutrino flavor transitions during the last decade. Even though it serves as the leading 
description, other mechanisms could be responsible for transitions on a subleading level. In 
this paper, we will develop an approximative framework for neutrino oscillations with so- 
called nonstandard interactions (NSIs) as such subleading effects. This framework is derived 
in the limit when the small mass-squared difference is negligible, i.e., Am^ — > 0, as well 
as when the effective matter interaction part of the Hamiltonian has only one eigenvalue 
significantly different from zero. The main effects of NSIs will be parametrized by v and f3, 
which will be defined in Sec. IIHI 

Recently, NSIs with matter have attracted a lot of attention in the literature. In some 
sense, the potential for these interactions can be viewed as a generalization of the coherent 
forward-scattering potential, which describes the effects of ordinary (or standard) interac- 
tions with matter 1L and indeed, gives rise to the famous Mikheyev-Smirnov-Wolfenstein 
(MSW) effect [l[ [2, 3]. Especially interesting are approximative frameworks that consider 
neutrino oscillations with NSIs which could be used for atmospheric neutrino data in the 

fin 

v^-fr sector, since these data are sensitive to interactions of tau neutrinos |4|, [5[. However, 
there are also other experimental situations, where one could investigate tau neutrino inter- 
actions such as with the MINOS experiment . In addition, in Ref. |9|, the authors 
have studied an approximative two-flavor neutrino scenario that is similar to our framework, 
but instead used for the OPERA experiment. Nevertheless, one should be aware of the fact 
that NSI effects are sublet effects and nave been exp—v contained Q Q 

The paper is organized as follows. In Sec. [Til we present the general scheme for three-flavor 
neutrino oscillations including NSIs. Next, in Sec. IIHI we develop an approximative two- 
flavor neutrino oscillation framework including subleading NSI effects as well as we compare 
our framework with earlier results. Then, in Sec. IIVI we analyze our framework numerically. 
In Sec. |Vl we study applications to experiments. Finally, in Sec. I VI I we summarize our 
results and present our conclusions. 
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II. NONSTANDARD INTERACTIONS IN NEUTRINO OSCILLATIONS 



Standard three-flavor neutrino oscillations in vacuum can be described by the following 
Hamiltonian in flavor basis: 

1 



#o = ^t/diag(0,Am^A m ^)f/t, 



2E 

where E is the neutrino energy, Aro^ = mf — is the mass-squared difference between the 
ith and jth mass eigenstate, and U is the leptonic mixing matrix. By adding the effective 
matter Hamiltonian 



(1) 



matter 



diag(v / 2G F A^ e , 0, 0) = Vdiag(l, 0, 0) 



(2) 



to the vacuum Hamiltonian Hq, standard matter effects on neutrino oscillations can be stud- 
ied. In order to investigate the effects of NSIs between neutrinos and other fermions, another 
more general effective interaction Hamiltonian can be added. This effective Hamiltonian will 
be of the form 



/ 
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NSI 
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(3) 



where 



-a/3 



Z^ fc «/3 AT ' 



Nf is the number of fermions of type /, and we have assumed an unpolarized medium. Thus, 
the full Hamiltonian is given by 



H — Hq + i^matter + -f^NSI, 

which describes three-flavor neutrino oscillations in matter including NSIs. 



(4) 



III. REDUCTION TO TWO-FLAVOR EVOLUTION 



Since the experimental bounds on e m and e ail are relatively stringent [lOj, we focus on 
the case where the interaction part of the Hamiltonian takes the form 



#int — -^matter + -f^NSI — V 



1 + See 
















e* 








yC/ NS idiag(i;,0,O^si. ( 5 ) 



where the unitary matrix 

' c p sp& 
Z/nsi = 10 

v -spe-* Cp J 

defines the matter interaction eigenstates, Cp = cos(/3), sp = sin(/3), and V is the standard 
MSW potential. In this framework, the NSIs are parametrized by the two matter interaction 
eigenvalues v and £ as well as the matter interaction basis parameters (3 and (b and we have 



£er = SpCpe^^ - v), e TT = vsl + ^ 



(6) 



Thus, the standard neutrino oscillation framework is recovered with (3 = £ = and v = 1. 

From the results of atmospheric neutrino experiments, we know that even high-energy 
Ufj, oscillate, leading to £ <C 1 0, 0] for the matter composition of the Earth. We will 
here consider the limit Vv ~ Amy (222) ^> V£, which is the applicable limit for standard 
neutrino oscillations with energies of a few GeV to tens of GeV inside the Earth, while 
earlier papers % |5] have studied the limit Vv ^> Am^/ (222) ~ V£. In the limit of £ — > 
and Am| 1 /(222) — > 0, the full three-flavor Hamiltonian takes the form (up to an irrelevant 
addition proportional to unity) 



22 = Aa/^a^ + Vva v al, 



(7) 



where 



f s 13 e- iS ^ 



S23 c 13 
C23 c 13 



and a v 



J 



C 





\ 



V spe-^ J 



are the third column of the leptonic mixing matrix U and first column of the matter inter- 
action mixing matrix U^si, respectively. Since this Hamiltonian is defined using only two 
linearly independent vectors, a a and a v , there must be a third linearly independent vector 
do for which Ha = and which is orthogonal to a a and a v . The components of this vector 
are given by 

a = ^J a A x cl v , (8) 

where d = cos(#') is a normalization factor with sm(9 r ) - - 
the product 



a A a v \ and we have introduced 



(A x B)i = s ijk A*B* k , 



(9) 



which is orthogonal to both A and B as well as antilinear in both arguments. With our 
parametrization, we have 



s 12 = (S13C/3 - c 13 Sf3C 2 3) 2 + c 2 3 sin(2/3) sin(26>i 3 ) sin 2 ' 



(10) 



Since ao is an eigenvector of the full Hamiltonian regardless of the matter density, its evo- 
lution decouples from that of the other two states and it will only receive a phase factor. 
Thus, it will be possible to describe the evolution of the two remaining states in an effective 
two-flavor framework. Indeed, if we choose the basis 



a x = e lctx a v , a y = e iay a v x a 



(11) 



(note that a y is already normalized as a v and ao are orthonormal), with a x and a y chosen 
such that a)^a x and a} A a y are real, then the two-flavor Hamiltonian in this basis takes the 
form 



H 2 = Vv 



1 \ Is' 2 s'd 



+ A 







s'd c' 2 



(12) 



which is an ordinary two-flavor Hamiltonian in matter with the vacuum mixing angle 9', the 
mass-squared difference Am^, and the MSW potential Vv. In order to compute the neutrino 
flavor evolution in this basis, any method applicable to two-flavor neutrino oscillations can 
be used. 

While the above framework makes it easy to reduce the three-flavor evolution to a two- 
flavor evolution in the given basis, the final results have to be projected onto the flavor basis, 
where 



M M M 





W 



1 

v / 







(13) 



The two bases are related by a unitary matrix V, which is given by 

V ai = v\di. 

With this notation and the full evolution matrix given by 



(14) 



S = e 



_i$ 











Sxx 


Sxy 





0* 
^xy 


O* 
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(15) 



in the Oj basis, the general neutrino oscillation probability becomes 



P 



a/3 



+ dvM' + \v Py v:s)p,„ 



xy 



+2 Re 
+2 Re 



Vp V: o e-^^p7 x (V* x V ax e iA * + V* y V ay e~ iA *) 



+2 Re [P xx V Px V: x V; y V ay e- 2 ^ - P xy V Px V: y V; y V ax e~ 2 ^] 



~\~2 \/ P xx P xy Re [Vp x V; y {\V ay \ 2 - \V ax \ 2 )e 



2\ -i(A x +Aj,) 



+2^Pj%Re [V: x V ay {\V Px \ 2 - \V Py \ 2 )e-^-^] , 



(16) 



where S xx = \/P xx e" lKx and S xy = P xy e~ lAy (with this notation P xx is the two-flavor 
survival probability and P xy is the two-flavor transition probability). 

From the above, it becomes apparent that the final neutrino oscillation probabilities are 
quite complicated unless further assumptions are made. This is mainly due to the fact that 
it is necessary to project the evolution matrix onto the flavor basis. In the case of (3 — 0, 
we regain the result from ordinary neutrino oscillations where a v = u e and the u e survival 
probability takes a very simple form, i.e., P ee = P xx . For simplicity, we present the analytic 
results only for the case 5 = = 0, where the matrix elements V ai are taken to zeroth order 
in S13. This simplification will not lead to any large errors, since V is fixed and S13 is small. 
However, we keep S13 7^ in the definition of 6', since there can be resonance effects in the 
two-flavor evolution. The resulting neutrino oscillation probabilities are then given by 

c'oPee — s 23 S /3 + C /3 {^xx + ^C^S^C^ \^xy (^-y) ~ Pxx Sin 2 (A x )] 

+ W P xxPx y c' s p c 23 (2c f2 - 1) sin(A x ) sin(A 2/ )| 
+2s 2 23 slcl I [cos(A ;E - $) - 244 sin($) sin(A ;E )] 



r' 4 P 

r' 4 P 

c r e/j. 



r f4 P 



+ 2^/P xy c' s l 3C23 sin(Aj,) sin($)| , 

C /3 C 23 Pxx s 23 ^ V ' Pxx C f3 C 23 S 23 cos {-^x + < ^ > )i 

+ 2slc 2 23 [P xx - y/P^cos(A x + $)] 
+ 2d^spC23\fPl y cos(A ?/ + $) - a/^cos(A ;e - A y )J | , 



+ SgC Q P Q 



I cJ4 (l- v / ^) 2 + 4 v / i^sin 2 ^^^^ 

+ 2a/P^c 2 9 s /3 c 23 c / cos(A j/ - $) - ^/p^cos(A x + A y ) 



xy 



(17) 
(18) 

(19) 
(20) 
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where Cg 2 = 1 — s|c| 3 . The five remaining probabilities can be trivially deduced from the 
above results using unitarity conditions. Note that Eqs. f|T7]) - fl20|) will be illustrated below 
in Figs. [Ml 



A. Comparison and consistency with earlier results 

As mentioned above, in Refs. , the limit of v — > oo and ~ Am 2 il /(2E) is 

studied. They also study the limit of £ — » 0, which becomes similar to the approach taken 



in Ref. 121 ]. where the limit of large matter effects was studied (the difference being that the 
former includes NSIs and the latter includes the effects of nonzero #13 and Amf-J. In this 
limit, the eigenstate a v decouples, and since V^ x = 0, the oscillates in a pure vacuum two- 
flavor system. It is found that the effective neutrino oscillation parameters in this system 
are given by 



Ami = Am^[c 2 ,(l + c 2 ) - S 2 ] 2 /4 + S 2 e c 2 , (21) 
tan(2ft m ) = 2 ^ 2 , (22) 

where C29 = cos(2# 23 ) and S29 = sin(26 l 23). 

In our framework, it is easy to obtain the corresponding quantities as 

Am 2 m = 2E{a [ y Ha y — ajifao) 

= Ami/ 2 , (23) 

tan(0 m ) = |V^ y |/|V^o| 
C13S23 



(24) 

C13C/3C23 + S/3S13 

where the last equality holds only when = 5 = and becomes more complicated when 
this is not the case. It is easy to show that these results are the same as those of Refs. 
in the limit of #13 — > 0, where we obtain 

Am 2 m = A^3i(l-4 C 2 3 ) ? (25) 
tan(U = (26) 

C/3 



which are slightly simpler forms of Eqs. (|21j) and ( 1221) . Note that Eqs. (1231) and ( 1241) are 
valid for any value of #13, whereas Eqs. ( 1211) and (1221) are only valid in the limit #13 — > 0. 

4 < Am^. 



Our results also agree with those of Ref. 



12 



in the limit of no NSIs and Am™ <C Am? 



Parameter Value 



Am§! 2.7 • 10" 3 eV 2 
Am^ 8 • 1(T 5 eV 2 



n-2 

h 3 



33.2° 
45° 
or 8° 



TABLE I: The neutrino oscillation parameters used in our numerical examples. The value of #13 
is specified in each figure. In addition, for all examples, we have used 5 = 0. 

TV. NUMERIC ANALYSIS OF OSCILLATION PROBABILITIES 

In this section, we present numerical results to show the accuracy of the two-flavor ap- 
proximation presented in the previous section. In order to perform this, we will present 
the neutrino oscillation probabilities in terms of neutrino oscillograms of the Earth, i.e., 
probability iso-contours in the zenith-angle-neutrino-energy plane. Since we are mainly in- 
terested in the accuracy of the approximation (rather than making quantitative predictions), 
we use a simplified Earth model, where the mantle and core are approximated as having con- 
stant matter densities (4.65 g/cm 3 and 10.2 g/cm 3 , respectively). For the standard neutrino 
oscillation parameters, we have used the values given in Table HI 

As will be observed in Figs. [THSl the approximations of Eqs. (TTTj) - (l2"0"j) are surprisingly 
accurate in describing the shapes of the oscillograms, the main source of error coming from 
the negligence of the small mass-squared difference Am^ (putting Am 2 x = in the numerical 
simulations would lead to differences between the approximations and the numerical results 
which would be hardly noticeable). 

We will also observe that, in all cases, the NSI terms dominate the high-energy behavior of 
the neutrino oscillation probabilities, as expected. However, in addition, it will be apparent 
that the effective two-flavor mixing angle 6' plays a large role in determining the qualitative 
behavior at energies around the resonance energy. This includes discussions on where the 
resonance appears (both in energy and baseline), as well as on other resonant effects such 



as parametric resonance 
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15 



16 



17 



181 ] . For energies lower than the ones displayed 



in Figs. (THSl our approximation will become worse due to the fact that the effects of Am^ 
will no longer be negligible. 



P = -ti/8 P = P = 71/8 




m m o— i m o 

o o o 

I I I 

cos(a) 



FIG. 1: The muon neutrino survival probability as a function of zenith angle and energy for 
different values of v and (3. The colored regions correspond to the numerical results (darker regions 
for higher probability), while the solid black curves correspond to the two-flavor approximation 
introduced in Sec. IIIII In this figure, the leptonic mixing angle #13 is set to zero. 



A. The muon neutrino survival probability Pfj,fj, 

In the next section, we will study the muon neutrino disappearance channel at a neutrino 
factory using the General Long Baseline Experiment Simulator (GLoBES). Thus, our first 
numerical example, shown in Fig. dj is the muon neutrino survival probability P MM . As was 
discussed earlier, at high neutrino energies, the muon neutrino will oscillate in a pure two- 
flavor vacuum scenario with the effective parameters given by Eqs. (1251) and ([261 . since we 
have used #13 = in the construction of Fig. [IJ Furthermore, the choice of #13 = and the 
approximation that Am^ is small imply that the low-energy neutrinos also oscillate in a 
two- flavor vacuum scenario, but with the parameters Am 2 = Am^ and 6 = 623- In the case 
of (5 = 0, these two cases are obviously equivalent and the matter interaction effects do not 



9 



change the situation at any energy (since the v e -v e element is the only nonzero element of the 
interaction contribution to the Hamiltonian) . However, for (3 7^ 0, there will be resonance 
effects in the region where A ~ Vv . As can be seen from Fig. [TJ these resonance effects are 
present mainly around E ~ 10 3 7 MeV for v — 1 and at energies appropriately displaced 
from this for v = 0.5 and v = 2 (i.e., E ~ 10 4 MeV and E ~ 10 3,4 MeV, respectively). Note 
that the colored regions in Figs. [THS] represent the "exact" (or numerical) three-flavor results, 
whereas the black solid curves are the approximate two- flavor results given by Eqs. (}TTj) - (j20p . 
In general, these two-flavor results are in excellent agreement with the exact three-flavor 
results, and therefore, they can be used as good approximations in analyses of neutrino 
data. 



B. The electron neutrino survival probability P e 



As discussed in Ref. [19], the electron neutrino survival probability P ee can be a very 
useful tool when considering ways of determining the small leptonic mixing matrix element 
Ue3- However, as has been discussed earlier in this paper (and also in previous papers 
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20 



21]), the presence of NSIs can significantly alter the interpretation of a positive 
oscillation signal in this channel. In particular, the effective two-flavor mixing angle 9' is 
given by Eq. ffTUj) and the effective matrix element f/ e3 by fy] 

2EV , , 

U e3 ~ U e3 + 

&er A 9 C23. \Z I ) 

Am 31 

In Figs. [2] and El we present the oscillograms for P ee using different values of (3 and with 
6*i3 = and 8°, respectively. It should be noted that s' ~ in the case when #13 = 8° 
and j3 = tt/16 (this is of course the reason for choosing #13 = 8°). As can be observed 
in these two figures, the high-energy (E > 10 4,5 MeV) behavior of the survival probability 
P ee is practically unaffected by the change in t9 13 . This is clearly to be expected, since 
the vacuum neutrino oscillation terms are suppressed by the neutrino energy E, while the 
matter interaction terms remain constant. Simply neglecting the vacuum oscillation terms, 
the expectation at high energy is 

l - P ee = sin 2 (2/3) sin 2 (VvL) , (28) 

which can also be obtained as a limiting case of Eq. ( TT7l) . 
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P = -7C/8 p = -71/16 




P = 7C/8 p = 371/16 




-1 -0.5 0-1 -0.5 

cos(oc) 



FIG. 2: The electron neutrino survival probability P ee as a function of zenith angle and energy for 
#13 = 0j v = 1, and different values of f3. The colored regions and solid black curves are constructed 
as in Fig. [TJ 

However, for energies in the resonance regime, the qualitative behavior of P ee is clearly 
dominated by the value of the effective two-flavor mixing angle 9'. This is particularly 
apparent in the mid-row panels of Figs. [2] and [31 While the oscillation probability in the 
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P = -Tt/8 p = -71/16 




P = 7C/8 p = 371/16 




cos(oc) 

FIG. 3: The same as Fig. El but for 13 = 8°. 

resonance region is kept small in both the mid- left panel of Fig. [2] and the mid-right panel of 
Fig. [3] (corresponding to 9' ~ 0), the resonant behavior at E ~ 10 3 ' 8 MeV in the mid-right 
panel of Fig. [2] and the mid-left panel of Fig. [3] (both corresponding to 8' ~ 8°) are also very 
similar. 
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C. The neutrino oscillation probability P t 



The oscillation of v e into v^, at neutrino factories has been thoroughly discussed in several 



papers (see Ref. [22j and references therein) and is commonly known as the "golden" channel. 
Thus, in Fig. HI we present oscillograms for the neutrino oscillation probability P efJi for 
different values of #13 and f3. Comparing this figure with Fig. [3j we can observe that, for 
#13 = 8° and (3 = 0, half of the oscillating u e oscillates into v^. This is a direct result of 
the leptonic mixing angle #23 being maximal and has been known for a long time (see, e.g., 



Ref. [23]). However, this is no longer true when f3 7^ 0. For the two top panels and the 
lower-left panel, there are regions where P efJL > 0.5; thus indicating that more than half of the 
u e oscillates into in these regions. It should be noted that there is an important caveat, if 
#23 deviates from maximal mixing, the linear combination of and v T , which v e oscillates 
into in the standard framework with nonzero #13, is not an equal mixture. Therefore, even 
the standard neutrino oscillation framework can allow for P efi > 0.5. However, what is 
important to note is that the ratio r = P efl /P eT is constant in the standard framework, but 
depends on both baseline and energy when (3 7^ 0. For example, if we consider neutrino 
oscillations with 9 13 = 8° and (3 = — 71/ 16 at cos(a) ~ 0.3, then r = at high energies 
due to the decoupling of the when the matter interaction terms dominate, while r > 1 
in the relatively low-energy region around the resonance. Hence, measuring the energy and 
baseline dependence of r could be a way of determining whether NSIs are present as long as 
the small mass-squared difference Ara^ is negligible. As in the case of the neutrino survival 
probability P ee , we can observe that there are practically no oscillations in the resonance 
region for the cases where 9' ~ 0. 



D. The neutrino oscillation probability P^ T 

In order to accommodate all of Eqs. (|T7l) -( !20l) in our numerical study, we also show the 
results for the neutrino oscillation probability P^ T in Fig. [5j Again, we can observe the fact 
that the oscillations for given values of 9' turn out to be very similar. Especially, the panels 
where 9' ~ are almost identical and very similar to pure vacuum neutrino oscillations. 
In particular, the effects of going from mantle-only trajectories to core-crossing trajectories 
at about cos(a) ~ 0.84 is not apparent in these figures, the reason being the fact that the 
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FIG. 4: The neutrino oscillation probability P e ^ as a function of zenith angle and energy for 
different values of #13 and (3. For this figure, we have used v = 1 and the colored regions and solid 
black curves are constructed as in Fig. [TJ 



Hamiltonian is diagonalized by the matrix (V a i) given in Eq. ( fT4l) for all values of VE. 

It should also be noted that P^ e = P e/U , since we use a symmetric matter density profile 
and no CP-violating phases in this section. Thus, the differences between the figures for the 
P m probabilities and those for the P^ T probabilities (except for the interchange of dark and 
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cos(a) 

FIG. 5: The neutrino oscillation probability P^ T as a function of zenith angle and energy for 
different values of #13 and (3. For this figure, we have used v = 1 and the colored regions and solid 
black curves are constructed as in Fig. [TJ 

light areas) are summarized in the figures for the P efl probabilities as P^ e + P^ + P Mr = 1 . 
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V. APPLICATION TO A FICTIVE NEUTRINO FACTORY 



In this section, we study the impact of NSIs on the sensitivity contours of a Active neu- 
trino factory experiment. In order to achieve this, we have used the General Long-Baseline 
Experiment Simulator (GLoBES) 24], |25j] and a modified version of the standard neutrino 
factory experimental setup included in the GLoBES distribution (essentially "NuFact-II" 



from Ref. 



26]), modified to only take the disappearance channel into account and to 



have a baseline of L = 7000 km (from Fig. [TJ we note that the NSI effects should be rela- 



tively strong at this baseline, it is also very close to the "magic" baseline [27J], which means 
that the impact of the solar parameters is expected to be small). This neutrino factory 
setup corresponds to a parent muon energy of 50 GeV and an assumed target power of 
4 MW. The running time is set to four years in each polarity and the detector is assumed 
to be a magnetized iron calorimeter with a fiducial mass of 50 kton. For the neutrino cross- 
sections, we have also used the files included in the GLoBES distribution, which are based 
on Refs. 28|, [29]. We consider neutrino energies in the range 4-50 GeV, divided into 20 
equally spaced energy bins. The simulated neutrino oscillation parameters are the same as 
in the previous section, i.e., they are given in Tabled! The simulated 6*i 3 has been put to 
zero for simplicity. 



A. Determination of standard parameters 



One important question when considering NSIs is how the inclusion of NSIs can alter 
the sensitivity of an experiment. This subject has been studied for a number of different 



inis suDiect 

QflflflQ 



neutrino oscillation experiments [4J, [5|, |6|, |7|, |9|, |20|, [21|, [30|, [31 ]. At the neutrino factory 
in question, the z/ M disappearance channel will normally be a very sensitive probe to the 
neutrino oscillation parameters Amjj and 823. However, with the inclusion of NSIs, there 
is an additional degeneracy, leading to experiments becoming less sensitive to the standard 
neutrino oscillation parameters. Note that this degeneracy has been investigated before in 



the case of atmospheric neutrinos in Refs. 



1. 



5|. In particular, the major part of the neutrino 



energies at the neutrino factory described before are well above the resonance energy. Thus, 
the induced degeneracy in the sin 2 (6 ) 2 3)-Am| 1 plane should be well described by Eqs. ( 1231) 
and (El. 
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sin 2 (0 ) 



FIG. 6: The sensitivity contours (2 d.o.f.) in the sin 2 (2^23 )~Am| 1 plane for the disappearance 
channel at the fictive neutrino factory experiment described in the text. The black curves corre- 
spond to the sensitivity if it is assumed that there are no NSIs present, while the colored regions 
correspond to the sensitivity if NSIs within current bounds are taken into account. The gray circle 
corresponds to the simulated parameters and the gray curve to the expected direction of the NSI 
degeneracy. The confidence levels are 90 %, 95 %, and 99 %, respectively. We have assumed 
Gaussian priors of e ee < 2.6, e eT < 1.9, and e TT < 1.9 at a 90 % confidence level. 



In Fig. [HI we can observe how this degeneracy manifests itself in the sensitivity of the 
disappearance channel. While the sensitivity contours are symmetric around sin 2 (6 l 23) = 0.5 
when not including NSIs (i.e., (3 = 0), they extend to smaller values of 9 2 3 and larger values 
of Arrigj when NSIs are included (i.e., (3 ^ 0). The conclusion that the sensitivity contours 
will extend to smaller #23 and larger Amf 1 is general and apparent from the form of Eqs. fl23|) 
and (I24p . where the same 8 m and AmJ, can be achieved for a class of different #23 and Amf^ 
given by 



tan(# 23 ) = C0taa(0 m ), (29) 

„2 2 
b /3 C 23 



Am 3i = T^fer; (30) 
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V 



FIG. 7: The sensitivity (2 d.o.f.) to the NSI parameters v and (3 for the disappearance chan- 
nel at the Active neutrino factory experiment described in the text. The dashed curves are the 
systematics-only sensitivities and the solid curves include parameter correlations with the standard 
neutrino oscillation parameters (with errors comparable to the current errors). The contours were 
constructed assuming no NSIs and the + corresponds to this situation. 



in the case when # 13 = 0. This degeneracy has also been marked in Fig. [6] for # 2 3 and Amf 
equal to the simulated values when j3 = 0. 



B. Determination of NSI parameters 

Another interesting question in relation to NSIs and future experiments is what the ex- 
perimental sensitivity to NSIs is, e.g., what bounds the experiments could put on the NSI 
parameters. Thus, in Fig. [7J we show the sensitivity contours in the v-f3 plane (the other 
NSI parameters have been marginalized) corresponding to the z/ M disappearance channel at 
a neutrino factory described above. From this figure, we notice that the given experiment is 
not very sensitive to the NSI parameters as long as the knowledge on the standard neutrino 
oscillation parameters has errors which are comparable to those of the current experimen- 
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tal bounds. Essentially, for large v (such that the entire energy spectrum is well above 
resonance), the bound on (3 comes from the mismatch between the effective 9 m and Am^ 
measured by the neutrino factory and the external measurements on # 2 3 and Am^ and the 
errors of both are propagated to the error in the determination of j3. However, if the external 
errors on the standard neutrino oscillation probabilities are very small (i.e., such that we 
do not need to consider correlations with the standard parameters) , then the bounds on the 
NSI parameters are significantly improved. It should be noted that the insensitivity to (3 for 
small \v\ is clearly due to the fact that it is not possible to determine (3 when the interaction 
effects are suppressed. We should also note that the given experiment is clearly not the 
best to use in order to distinguish NSIs from standard neutrino oscillations. For example, 
it should be more viable to detect NSIs through other oscillation channels or through the 
combination of the disappearance channel with one of the appearance channels. 

VI. SUMMARY AND CONCLUSIONS 

We have considered the effects of NSIs on the propagation of > GeV neutrinos in the Earth 
with the assumption that Am^ is small and that e ee , e eT , and e TT are the only important 
NSI parameters. We have shown that, if the NSI parameters are constrained to be on the 
"atmospheric parabola," then the neutrino flavor evolution matrix can be straightforwardly 
computed using a simple two-flavor scenario in matter with a mixing angle given by Eq. (fTOl) . 
This two-flavor scenario is the natural NSI extension of the usual two-flavor approximation 
used in standard neutrino oscillations when Am^ is small. While this two-flavor scenario, 
unlike the standard two-flavor scenario, is not diagonal in flavor space, it should be noted 
that if NSI effects are present in the propagation, then they should also be present in the 
creation and detection processes, meaning that the actual flavor space is no longer equivalent 
to the weak interaction flavor space. Deviations from the atmospheric parabola or other NSI 
parameters can be easily treated as a perturbation to the two-flavor scenario. 

The impact of NSIs on neutrino oscillations in the Earth as well as the accuracy of our 
two-flavor approximation have been studied using neutrino oscillograms of the Earth. It was 
noted that the resonance regions in the neutrino survival channels are mainly dependent on 
the mixing angle of the effective two-flavor scenario, implying a degeneracy between # 13 and 
the NSI angle (3. This degeneracy is preferentially broken by using data from above or below 
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the resonance, or by studying the flavor combination into which neutrinos actually oscillate. 

Finally, we have given an example of the impact of the NSIs in a Active neutrino fac- 
tory with a baseline of L = 7000 km studying the muon neutrino disappearance channel. 
We have shown that the degeneracy between the standard and NSI parameters extends in 
the expected direction and computed the actual sensitivity to the NSI parameters, both 
by assuming no degeneracy with the standard parameters (i.e., the standard parameters 
assumed to be exactly known) and by assuming correlations with the standard parameters 
(in this case, Gaussian priors with a size corresponding to the current experimental errors 
were used). Again, it should be stressed that this part is only intended to show the impact 
of the NSI parameters in a Active experiment and does not offer a realistic approach for 
searching for NSIs. Such considerations can be found in Ref. jsoj]. 
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